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Abstract 

Based on an extended compound nucleus model, isospin effects in statistical frag- 
ment emission from excited nuclear systems are investigated. An experimentally 
observed scaling behavior of the ratio of isotope yields Yi(N,Z) from two similar 
emitting sources with different neutron-to-proton ratios is predicted theoretically, 
i.e., the relationship of Y%jY\ oc exp(aN + (3Z) is demonstrated. The symmetry 
energy coefficient C sym extracted from the simulation results is ~ 27 MeV which 
is consistent with realistic theoretical estimates and recent experimental data. The 
influence of the surface entropy on the isoscaling behavior is discussed in detail. It is 
found that although the surface entropy increases the numercial values of isoscaling 
parameters a and 0, it does not affect the isoscaling behavior qualitatively and has 
only a minor effect on the extracted symmetry energy coefficient. 
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1 INTRODUCTION 



The search for a nuclear equation of state (EOS) has been one of the dominant 
factors driving interest in heavy-ion collisions at intermediate and high ener- 
gies, ever since such beams have become available. Exploration of the isospin 
dependence of the EOS is an inseparable part of the task. Originally, the 
latter studies were confined to use stable projectile beams with only a mod- 
erately wide range of isotopes between neutron-poor and neutron-rich nuclei. 
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However, exotic secondary beams of projectile nuclei with extreme neutron- 
to-proton ratios that have become available recently offer new and intriguing 
opportunities to study isospin physics in heavy-ion collisions [Tf2"f3"] . These 
studies have already resulted in the discovery of several interesting isospin- 
related systematics [4.5.6.7J, prompting further efforts both, in experiments 
and in theoretical modeling. In particular, a so-called isoscaling analysis of 
light fragment emission in the decay of very hot nuclear systems produced in 
heavy-ion collisions has attracted attention as a possible tool for deducing the 
nuclear symmetry energy from the relative fragment yields [7 , 8 , QfTTVTU] . 

For strongly damped collisions fragment isotopic yields were found to follow 
a u Q gg systematics" [T2lll3fl4j . Isotopic scaling, termed "isoscaling" [7118119] . 
is observed for other types of reactions such as evaporation [71fl~5] . fission 
[HUT], projectile fragmentation [THUS], and multifragmentation [7115115] . This 
scaling law refers to a general exponential relation between the ratios of yields 
Yi(N, Z) for fragments (N, Z) emitted from systems which differ only in their 
isospin or {% L )i- In particular, if two reactions lead to primary systems i — 1 
and 2 having approximately the same temperature but different isospin, the 
ratio B,2i(N,Z) of the experimental yields of a given fragment (N,Z) emit- 
ted from these systems exhibits an exponential dependence on the fragment 
neutron number N and atomic number Z of the form, 



Here, a and (3 are the isoscaling parameters and C is an overall normalization 
constant. 

On the theoretical side, isoscaling has been extensively examined in the anti- 
symmetrical molecular dynamics model [20] , a Boltzmann-Uehling-Uhlenbeck 
model [5], a lattice gas model [2"Tf2"2"] . the expanding evaporating source model 
[23] and in statistical multifragmentation models [8ll9llllf24] . In the present 
paper it will be shown how the essential features of isoscaling observed in 
experiment are related to the nuclear symmetry energy in an extended com- 
pound nucleus (ECN) model. This model is closely related to that known 
from fission studies, and its central notion is a relatively high entropy (par- 
ticularly at moderately high excitation energies) associated with the diffuse 
nuclear surface region (as opposed to bulk matter). Application of this ECN 
model reveals a new mechanism of nuclear fragmentation caused by the soft- 
ening of the diffuse nuclear surface at high excitations [25J. The ECN model is 
equivalent to conventional compound nucleus models at low excitation energy 
as represented, e.g., by GEMINI [26], but greatly extends the validity of the 
compound nucleus concept towards high excitations [27]. Therefore, the ECN 
model provides a unified description of statistical emission of light particles 
and fragments in a wide range of excitation energies. While the present model 
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is still somewhat schematic, it is based on a physically transparent picture and 
thus provides direct insight into the phenomena of interest. 



2 THEORETICAL FRAMEWORK 



The probability p of emitting a fragment from an equilibrated compound nu- 
cleus (CN) is evaluated using the Weisskopf formalism 



p (X e AS = e Saaddls ~ Seq , (2) 



where S eq and S sa ddie are the entropies for the equilibrated CN and a saddle- 
point configuration of touching spheres, respectively. Within the Fermi gas 
model, the entropies for the two configurations can be calculated as 

S eq = 2\Ja A E^ ot , (3) 



and 

Ssaddie = S res + S frag = 2v/ ((Ires + a frag)E* a d dle . (4) 



In Eqs. (3) and (4), a a, a res and aj rag are the level density parameters of the 
CN system at equilibrium, of the residue, and of the fragment, respectively. 
The part Sg of the entropy S of the system, associated with the diffuse surface 
domain, has been found [29] to have a pronounced effect on the fragment 
emission probability. One has 

S = Sy + S s , (5) 



where Sy is the entropy of bulk matter. Consequently, the level density pa- 
rameter a includes volume and surface terms 1 30 1 , 



a v + a s = a v A + a s A 2/3 F 2 , (6) 



Here A is the atomic number, F 2 is the surface area relative to a spherical 
shape. E*aj dle in Eq. (4) is the thermal excitation energy of the system in the 
saddle-point configuration. The latter quantity is calculated as 

Efaddle — Etot ~ ^saddle, (7) 
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where V sa ddie is the collective saddle-point energy. 



For the level density parameter a, the parametrization of Toke and Swiatecki 
[30] was employed with a v = 1/14.6 MeV" 1 and a s = 4/14.6 MeV" 1 . The 
calculations assume saddle-point shapes to be represented by two touching 
spheres. V sa ddie is the difference in deformation energies for the saddle-point 
shape and equilibrium-state shape. It contains contributions from Coulomb, 
volume, and surface energies. Volume and surface energies depend on the 
system isospin / (defined as (N — Z)/A) in a functional form —a v (l — n v I 2 )A 
and a s (l — k s I 2 )A 2 ^ 3 , respectively. The parameters are taken from [21] as 
a s = 21.13, k s = 2.3, a v = 15.9937, k„ = 1.927. The nuclear temperature can 
be obtained from the commonly used Fermi-gas model relationship between 
the temperature T and the excitation energy of the system E^ ot : 



The present calculations do not account for the effects of an expansion of the 
CN prior to its decay, so the nuclear density is that of the ground state. There- 
fore the additional influence of the nuclear matter density on the extraction of 
symmetry energy coefficient is not studied here. The present work concentrates 
on the temperature dependence of the relation between symmetry energy and 
isoscaling parameters. 



3 RESULTS AND DISCUSSION 



In this work several pairs of equilibrated CN sources with proton number 
Z s = 75 and mass numbers A s = 165, 175, 185, and 195, are considered at 
initial excitation energies of E^ ot /A = 2, 3, 4, 5, 6, and 7 MeV/nucleon. As the 
ratio B,2i(N, Z) is insensitive to sequential decay, it carries information on the 
original excited fragments prior to their decay [18]. In particular, it has been 
found that the values of a and f3 are not much affected by sequential decay of 
the primary fragments [9]. These observations justify the neglect of sequential 
decay in the following analysis. 

The yield ratio B,2i(N,Z) is constructed using the convention that index 2 
refers to the neutron- rich system and index 1 to the neutron-poor one. The 
observation summarized in Eq. (1) shows that experimental yield values of 
ln(R 2 i(N, Z)) plotted vs. N [Z fixed) or Z (N fixed), produce straight lines. 
This feature also emerges from the present MCN model. As a demonstration, 
Fig. 1 shows the theoretical yield ratios ln(R 2 i{N, Z)) plotted vs. fragment 
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neutron number N, for individually values of Z = 6 — 9 (top panel), and vs. 
proton number Z, for individually values of N = 6 — 9 (bottom panel). The 
ratios have been calculated for the two source pairs (A 2 = 175, Ax = 165) and 
(A 2 = 185, A\ = 165) at an excitation energy of E^ ot /A = 5 MeV/nucleon. 
One can see from Fig.l that isotope ratios of the same element Z, or isotone 
ratios of the same N, tend to lie on a logarithmic straight line. The solid lines 
represent linear fits to each series of Z isotopes and each series of N isotones, 
respectively. It is obvious that these fit lines are nearly parallel to each other. 
The positive slope in the upper panel of Fig.l illustrates that neutron-rich 
fragments are more easily produced from the more neutron-rich sources, as 
can be expected. Analogously, the negative slope in the bottom panel of Fig.l 
illustrates that proton-rich fragments are more readily produced from the more 
proton-rich emitting sources. These are the typical features associated with 
"isoscaling" behavior showing that the scaling behavior of the isotope and 
isotone yield ratios appears naturally within the present CN model framework. 

It was reported [29] that the ratio of a certain fragment (N, Z) emission 
probability from two equilibrated systems can approximately be written as 
P 2 (N,Z)/Px(N,Z) oc exp[(Vx(N,Z) - V 2 (N } Z))/T], where V(i = 1,2) is of 
the same meaning as V sa ddie, namely the interaction energy at the saddle point. 
An analysis indicates that in the difference (Vx — V 2 ), most terms cancel ex- 
cept for terms that are directly related to the isospin of the CN system. These 
isospin-related terms can be combined as (a v K v — a s n s / A 1 ^) (N — Z) 2 / A = 
C sym (N - Z) 2 /A = E sym . Here C sym (= a v n v - a s K s /A 1/3 ) is symmetry en- 
ergy coefficient |32j incorporating both the volume and surface contributions 
to the symmetry energy E sym [9]. Hence the value of (Vx — V 2 ) is dominantly 
determined by the difference in symmetry energy between the two emitting 
sources and their residues. This demonstrates that the origin of the isoscaling 
phenomenon found in the present frame can be traced to the symmetry energy 
term in the nuclear binding energy, as also suggested in other work (e.g., [7]). 

The isoscaling parameters a and (3 can be extracted from fits of model pre- 
dictions to the data points shown in Fig.l. An average value of a is calculated 
over the range 5 < Z < 9, an average (3 is calculated over the range 5 < N < 
9. Figure 2 depicts the dependence of the average coefficients a and (3 vs. 
excitation energy (left panel) and vs. inverse of the temperature (right panel). 
One notices that the absolute values of a and | (3 | decrease monotonically 
with excitation energy, implying that isospin effects decrease with increasing 
excitation energy. There is a significant sensitivity to the excitation energy 
at low energy, but both the sensitivity to excitation energy and the overall 
isospin effect are weakened at very high excitation. In addition, it is rather 
evident that both a and (3 show a linear dependence on 1/T. 

The surface contribution as to the level density parameter a appreciably en- 
hances the fragment emission probability via significant surface entropy effects 
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[29] . It is therefore interesting to examine the effects of surface entropy on the 
isoscaling phenomenon in the present model. Results are displayed in Figs. 3 
and 4, to be compared to Figs.l and 2 which have been calculated without 
surface contribution to the entropy. Evidently, familiar isoscaling behavior is 
predicted well in either case. The linear relationship between isoscaling pa- 
rameters and inverse temperature is not affected by the surface contribution 
either, as seen from the right panel of Fig.4. However, an account of the surface 
term predicts larger values of a and (3 (comparing Fig.4 and Fig. 2), implying 
that the surface entropy effect on fragment emission becomes stronger with 
increasing isospin of the emitting sources, also as expected. 

Since the scaling behavior of ratios of fragment isotopic yields measured in 
separate nuclear reactions has been utilized to probe the symmetry energy 
[71ll0f33f32j . in the following a symmetry energy coefficient is extracted from 
the theoretical isoscaling systematically analytical expressions for a and (3. 
It is worth noting that, for a given E^ ot /A these four sources have the same 
temperature when the surface contribution to the level density parameter is 
neglected. 

It has been shown [91TTT] that the isoscaling parameters a and (3 are related to 
the symmetry energy coefficient C sym as 

« = 4% 1 [(t)!-(t)11 ( 9 ) 



and 



n n 



,N, 
A 



S )l] 

s 



(10) 



Equations (9) and (10) have been used to constrain the symmetry energy co- 
efficient C sym based on experimental data (e.g. Ref.JTU]). Figure 5 depicts the 
product a ■ T as a function of (Z s /A s )\ — [Z s j A^)\ and (3 ■ T as a function of 
{N s /A s )\ - (N s /A s )l for the initial four source pairs at various excitation ener- 
gies. All these systems with different source sizes and isospin asymmetries lie 
along one single straight line, which illustrates that the isoscaling parameters 
a and (3 are not sensitive to the system size. By fitting the theoretical data in 
Fig. 5 with Eqs. (9) and (10), a symmetry energy coefficient C sym can be ob- 
tained to be 27.3 ± 0.1 MeV, which is close to the standard liquid drop model 
value C sym = 25 MeV and also in reasonable agreement with that obtained in 
Ref.[TS]. In that experimental work [18J, a value of C sym = (27.2 ± 2.2) MeV 
was deduced based on an experimental analysis of 25-MeV/nucleon 86 Kr + 
i24,ii2g n reac ti ons _ The present analysis demonstrates consistency with other 
realistic theoretical estimates [3413] and indicates the validity of an interpre- 
tation of isoscaling data in terms of the symmetry energy. 
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Although the inclusion of the surface term as in the level density parameter 
leads to a difference in the temperatures for the four systems, calculations 
indicate that for each excitation energy E% ot /A, the temperature differences 
among the four sources are less than 0.1 MeV, so the temperature value of four 
systems can be regarded as approximately constant when the excitation energy 
per nucleon is fixed. Figure 6 displays the effect of the surface entropy on the 
extracted symmetry energy coefficient. The value deduced from the theoretical 
data in Fig. 6 is C sym = (27.6 ± 0.1) MeV, which is comparable with that 
predicted without considering the surface entropy effects. The reason is that 
although the surface entropy increases the values of a and j3, it also decreases 
the value of temperature due to an increased level density parameter a [see 
Eq. (6)]. Consequently, a-T and (3-T are little affected by the surface entropy. 
This result indicates that surface entropy probably has a minor influence on 
the extraction of the symmetry energy coefficient. Moreover, it also implies 
that the isoscaling observable is a robust probe of the symmetry energy. 



4 SUMMARY AND CONCLUSIONS 



In conclusion, in the framework of an extended compound nucleus model, 
isoscaling behavior and its relation to the nuclear symmetry energy is revealed. 
The symmetry energy coefficient is found to be C sym ~ 27 MeV from an anal- 
ysis of theoretical "data", suggesting that isoscaling data can be interpreted 
in terms of the symmetry energy. In addition, in this work surface entropy 
influences on the isoscaling phenomenon have been studied. The present ECN 
model approach leaves sufficient room for further improved treatment of in- 
teresting physical effects [29], such as nuclear expansion. Such a more detailed 
analysis is required to study the evolution of the symmetry energy with the 
excitation energy, as has been explored recently also by Shetty et al. [55] . 
The model will be utilized to deduce the effects of nuclear expansion on the 
isoscaling parameters and will correspondingly illustrate how to deduce the 
density dependence of the symmetry energy in such more general and realistic 
scenarios. Work along this direction is in progress. 
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Fig. l. (Color online) The logarithm of the ratio of elements Z = 
5-9 isotopes yields (top panel) and of N = 6-9 isotones yields 
(bottom panel) from pairing source A 3 = 175 and A s = 165 (left 
column) as well as from pairing source A s = 185 and A s = 165 
(right column) at excitation energy E^ JA = 5 MeV/nucleon. 
Here the calculations do not consider the surface entropy, i.e., 
not including as in the expression for the level density parameter 
a [see Eq. (6)]. Solid lines are the linear fitting to the data points. 
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Fig. 2. (Color online) Dependence of isoscaling parameters a and 
| (3 | on excitation energy (left panel) and the inverse temper- 
ature (right panel) for various source pairs. Symbols in the fig- 
ures are a (solid symbols) or | (3 | (open symbols) from four 
source pairs 1^=175/^=165 (squares), Ya s =W1a s =165 (circles), 
1a=195/1a s =175 (up-triangles), and 1^=195/1^=185 (down-trian- 
gles). Here the calculations do not consider the surface entropy, 
i.e., not including as in the expression for the level density pa- 
rameter a [see Eq. (6)]. 
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Fig. 3. (Color online) Same as Fig.l but consider the surface en- 
tropy effect via inclusion of as in the expression for the level 
density parameter a [see Eq. (6)]. 
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Fig. 4. (Color online) Same as Fig. 2 but consider the surface en- 
tropy effect via inclusion of as in the expression for the level 
density parameter a [see Eq. (6)]. 
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A(N / A ) 2 




Fig. 5. (Color online) a-T (positive parts) and (3-T (negative parts) 

as a function of (Z s /A s )\ - (Z s /A s )\ or as a function of (N s /A s )\ - 

(N s /A s )\ of the sources for four source pairs, i.e., Ya s =i7s/Ya s =w5, 

y^=i85/yk,=i65, YA.=m/YA.=m, and 1^=195/^=185, at excita- 
tion energies of E^ ot /A = 2, 3, 4, 5, 6, and 7 MeV/nucleon. Here 
the calculations do not consider the surface entropy, i.e., not in- 
cluding as in the expression for the level density parameter a [see 

Eq. (6)]. 
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Fig. 6. (Color online) Same as Fig. 5 but consider the surface en- 
tropy effect via inclusion of as in the expression for the level 
density parameter a [see Eq. (6)]. 
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